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ABSTRACT 


The  main  objective  of  this  program  is  the  investigation  of  the 
pulsed  electrical  CO  laser  to  determine  the  relative  scaling  param- 
eters and  the  usefulness  of  high  power  CO  lasers  in  systems  applications. 
The  current  contract  has  included  measurements  of  small  signal  gain  on 
varying  vibrational  and  rotational  transitions,  operation  with  an  unstable 
resonator  and  the  design  fabrication  and  assembly  of  an  entirely  new 
transverse  flow  system  together  with  a new  heat  exchanger. 

Small  signal  gains  of  up  to  4 l/2  % per  cm  were  measured.  Maxi- 
mum gain  was  observed  on  the  9 P ( 1 0)  and  the  8P(10)  transitions.  The 
measurements  were  compared  with  parametric  performance  predictions 
derived  from  kinetic  modelling  calculations.  The  experimental  data 
agreed  qualitatively  with  the  modelling  calculations.  Some  measurements 
on  low  transitions  in  particular,  6 P(10)  and  5 P(13)  support  the  hypothesis 
that  resonant  absorption  by  overlapping  transitions  can  seriously  suppress 
gain  on  certain  lines.  This  mechanism  will  play  an  important  role  in 
determining  the  spectral  distribution  and  hence  the  atmospheric  trans- 
mission properties  of  the  high  pressure  carbon  monoxide  electric  discharg 
laser. 

The  new  flow  system  discharge  cavity  and  heat  exchanger  have 
been  designed  with  the  aim  of  obtaining  a near  diffraction  limited  beam 
from  the  pulsed  CO  laser.  The  design  considerations  are  described  in 
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I.  INTRODUCTION 


The  CO  laser  program  is  directed  toward  the  investigation  of  the 
performance  characteristics  of  the  pulsed  electrical  CO  laser  using  an 
electron  beam  sustained  discharge.  The  experimental  investigation  is 
aimed  at  evaluating  the  scaling  parameters  of  high  power  electrical  CO 
lasers  and  their  potential  systems  applications.  The  overal'  program 
includes  theoretical  modelling  to  provide  support  for  the  design  of  the 
experiment  and  the  interpretation  of  the  data. 

The  following  areas  of  investigation  were  included  in  the  current 
contract. 

1)  Measurements  of  the  temporal  variation  of  the  small  signal 
gain  over  the  entire  range  of  accessible  transitions. 

2)  Operation  with  an  unstable  resonator  over  the  central  region 
of  the  original  discharge  cavity. 

3)  Design  and  construction  of  an  intracavity  water  vapor  absorp- 
tion cell. 

4)  Design  ana  construction  of  a new  flow  system/dis char ge 
cavity  and  heat  exchanger  for  operation  over  the  temperature 
from  80°K  to  300°K  and  gas  densities  of  the  order  of  1 amagat 
with  the  design  goal  of  achieving  better  than  1 l/2  times  the 
diffraction  limited  beam. 

The  small  signal  gain  measurements  are  described  in  detail  in 
the  present  report  and  are  compared  with  parametric  performance 
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predictions  derived  from  the  AERL  kinetic  model.  Maximum  gains  of 
4 1/2%  per  cm  were  measured  on  two  transitions.  Anomalous  behavior 
was  observed  on  several  low  lying  transitions.  These  transitions  exhibited 
oscillating  gain/absorption  behavior  or  in  some  cases  only  absorption 
where  gain  had  been  predicted.  The  measurements  can  be  interpreted  when 
allowance  is  made  for  overlapping  transitions  which  results  in  gain  suppres- 
sion on  several  lines.  These  overlaps  will  play  an  important  role  in 
determining  spectral  distribution  of  the  high  pressure  carbon  monoxide 
electric  discharge  laser  and  will  influence  the  atmospheric  transmission 
properties  of  the  system.  Lasing  operation  with  the  unstable  resonator 
cavity  indicated  nonuniform  energy  deposition  across  the  discharge. 

The  new  flow  system/discharge  cavity  has  been  designed  for 
transverse  flow  operation  directed  vertically  upwards.  This  is  to  ensure 
that  buoyancy  effects  are  favorable  and  to  minimize  the  boundary  layer 
growth  on  the  room  temperature  side  walls.  A porous  plate  heat  exchanger 
system  has  been  designed  and  constructed  with  large  thermal  mass  to 
provide  the  order  of  10  sec  flow  operation  of  the  mass  flow  of  150  grams 
per  sec.  The  discharge  cavity  is  designed  to  be  compatible  with  the 
existing  20  cm  x 100  cm  electron  beam. 
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II.  GAIN  MEASUREMENTS 


Measurements  of  small  signal  gain  have  been  performed  in  the 

20  liter  E-beam  sustained  CO  laser  and  the  results  compared  with  para- 

13) 

metric  performance  predictions'  7 derived  from  kinetic  modelling  calcula- 
tions. ^ Maximum  gain  of  4 1/2%/cm  was  measured  on  the  9P(10)  and  8P(10) 
transitions.  The  measurements  support  the  hypothesis  ^ that  resonance 
absorption  by  overlapping  transitions  is  the  mechanism  responsible  for 
gain  suppression  on  certain  lines.  This  mechanism  will  play  an  important 
role  in  determining  the  spectral  distribution,  and  hence  the  atmospheric 
transmission  properties  of  the  high  pressure  carbon  monoxide  electric 
discharge  laser. 

A.  TECHNIQUE 

The  technique  is  conceptually  simple.  A C.  W.  CO  probe  laser 
beam  was  transmitted  through  the  active  discharge  volume  of  the  20  liter 
E-beam  sustained  CO  laser.  The  cavity  mirrors  were  removed  to  pre- 
vent  oscillation.  The  laser  was  operated  at  an  initial  gas  temperature 
of  1 30°K,  a density  of  l/3  Amagat  and  with  a gas  composition  20%  CO  = 

80%  N2.  The  electric  field  initially  13  KV/cm/atm  decayed  to  approxi- 
mately 10  KV/cm/atm  by  the  end  of  the  60  fxsec  current  pulse.  The 
voltage  droop  was  associated  with  discharging  the  capacitor  bank  supplying 
the  energy  to  the  electric  field. 
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B.  PROBE  LASER  CONSTRUCTION 


Considerable  attention  was  paid  to  the  design  of  the  probe  laser  to 
ensure  high  output  stabili*  T.  Details  of  the  probe  laser  construction  are 
shown  in  Figure  1.  The  laser  consisted  of  a 150  cm  long  thick  walled 
pyrex  tube  with  a 1.  25  cm  i.  d.  The  electrodes,  formed  from  high  purity 
nickel,  were  mounted  in  side  arms  and  located  one  near  each  end  and  one 
at  the  center  of  the  discharge  tube. 

Special  care  was  taken  in  mounting  the  electrodes  to  prevent  the 
discharge  seeing  the  edges  of  the  nickel  cylinders,  thus  reducing  the 
formation  of  hot  spots  and  significantly  improving  the  stability  of  the  dis- 
charge. This  design  also  suppresses  sputtering  from  the  outer  electrode 
surface  onto  the  glass  enclosure  reducing  chemical  degradation  of  the  gas 
composition.  The  flowing  gas  mixture,  He:N2:CO,  10:1:0.3  torr  was 
pumped  by  a 15  c.  f.  m.  fore  pump. 

The  conventional  method  of  providing  liquid  nitrogen  cooling  to  the 
discharge  tube  by  circulating  the  coolant  in  a concentric  glass  jacket  was 
not  adopted.  This  technique  may  degrade  the  performance  of  the  laser 
due  to  mechanical  vibration  associated  with  circulating  the  coolant. 

Temperature  variation  along  the  discharge  tube  may  also  occur  as  a con- 
sequence of  gas  pressure  gradients  in  the  jacket  introducing  variations  in 
the  temperature  of  the  liquid  nitrogen. 

Instead,  the  discharge  tube  was  immersed  in  a large  volume  bath 
of  freely  evaporating  liquid  nitrogen  which  obviated  both  problems  associ- 
ated with  the  concentric  cooling  jacket  scheme.  The  bath  consisted  of  a 
1/2  mil  thick  aluminized  mylar  foil  bag  sealed  onto  the  ends  of  the  tube  with 
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a low  temperature  epoxy  cement.  The  bag  was  formed  and  supported  by  a 
rectangular  wooden  box  insulated  with  a styrofoam  lining.  The  ends  of  the 
discharge  tuoe  exterior  to  the  bath  were  terminated  with  Brewster  angle 
calcium  fluoride  windows.  Vacuum  integrity  was  maintained  by  viton  0-ring 
seals  between  aluminum  supports  and  the  discharge  tube  and  windows. 
Heater  tapes  wound  around  the  discharge  tube  between  the  bath  and  the 
window  holders  prevented  excessive  cooling  of  the  0-ring  seals  and  main- 
tained the  windows  at  room  temperature  to  avoid  problems  of  condensation 
The  laser  cavity  was  formed  by  a 4 meter  radius  95%  transmitting  germa- 
nium mirror  and  an  original  grating  blazed  at  5.  4^  with  360  lines/inch. 

Both  elements  were  placed  in  gimbal  mounts  with  micrometer  tilt  adjust- 
ments and  equipped  with  piezo  electric  drives  to  provide  longitudinal 
mode  tuning.  The  mounts  were  each  bolted  to  massive  invar  plates  sepa- 
rated at  their  four  corners  by  2"  diameter  invar  rods,  thus  providing  a 
high  degree  of  thermal  and  mechanical  stability  for  the  laser  cavity.  The 
invar  frame  was  further  supported  on  a 12"  x 3"  steel  channel  which  also 
carried  ground  steel  rods  which  formed  on  optical  bench  for  supporting 
the  aluminum  window  mounts.  Finally,  the  channel  was  mechanically  iso- 
lated from  the  surroundings  with  an  auto- vibration  mounting. 

The  laser  was  excited  by  a dc  discharge.  Power  was  provided  from 
a current  regulated  power  supply  typically  operating  at  1 5 kV  and  supplying 
10  mA  through  500  kQ  ballast  resistors  connected  to  each  cathode. 

Single  transverse  mode  oscillation  was  achieved  using  an  intra 
cavity  6 mm  diameter  aperture  placed  close  to  the  diffraction  grating. 

Beyond  the  cavity  the  optical  path  was  defined  by  a 6 mm  diameter  iris  after 
which  typical  line  intensities  of  100  - 300  mW  were  obtained. 
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C.  SPECTRAL  OUTPUT 

lhc  wavelengths  of  the  lasing  transitions  were  measured  with  a 1/2 
meter  Jarrel-Ash  monochromat _,r  equipped  with  a 295  lines/mm  grating 
blazed  at  3.  7 p.  used  in  conjunction  with  a Au-Ge  Detector.  The  mono- 
chromator was  calibrated  with  mercury  lines  in  high  order.  Using  100  p 
slits,  the  accuracy  of  the  instrument  at  2000  cm"1  was  + 0.  5 cm"1. 

The  wavelengths  of  50  P branch  transitions  for  the  first  25  vibra- 
t-onal  levels  were  calculated  using  the  sxpression  for  .he  energy  levels  of 
the  anharmonic  oscillator  and  the  values  of  the  spectroscopic  constants  for 


the  CO  molecule  recently  determined  by  Eng  et  al. 


(4) 


The  line  separation  within  a vibrational  band  is  typically  4 cm"1, 

the  calculated  wavelengths  are  believed  correct  to  within  10'4  cm"1.  Thus, 

correct  line  identification  is  normally  assured.  Occasional  difficulties  are 

associated  with  overlapping  lines  encountered  at  the  band  extremities; 

however,  this  difficulty  could  be  resolved  by  noticing  variation  between 

the  spacings  of  lines  belonging  to  different  bands.  The  identification  was 

confirmed  by  predicting  the  identity  of  two  missing  transitions  in  the  8-7 

band  which  are  very  strongly  absorbed  by  overlapping  water  vapor  transi- 
(5) 

tions  and  thus  do  not  propagate  through  the  atmosphere.  Oscillation 

occurred  on  approximately  six  P branch  transitions  P(9)  - P(15)  within  the 
ten  vibrational  bands  V'-V  14^13  to  5^4.  The  calculated  frequencies  listed 
in  Table  1 correspond  to  those  transitions  which  were  identified  to  within 
the  precision  of  the  experimental  measurements  corrected  to  vacuum. 

The  small  signal  gain  measurements  were  made  by  monitoring  the 
probe  laser  intensity  with  a Au-Ge  detector  after  propagation  through  the 
discharge.  The  temporal  variation  of  the  detector  signal  during  and  after 
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TABLE  I 

SPECTRAL  OUTPUT  OF  THE  PROBE  LASER  AND  CALCULATED 
FREQUENCIES  OF  THE  LASING  TRANSITIONS 


Transition 

r 

. V p / T \ 


Calculated 

Frequency 


Trans  ition 
1 


1977  • 

■ 2738 

1973  « 

■ 2955 

1969  ■ 

2837 

1965  • 

2386 

1961  • 

16  02 

1957  • 

0486 

195!  • 

455  1 

1947  • 

51  1 9 

.943  • 

5352 

939  • 

5250 

935  • 

4817 

931  • 

4052 

927  • 

2958 

Calculated 

Frequency 


92  1 

8026 

913 

8858 

905  • 

8361 

896  • 

16  93 

892  • 

2626 

888  • 

3225 

884  . 

3493 

880  . 

3429 

876  • 

3036 

1870  • 

6134 

866  • 

7417 

86  2 • 

8367 

858  • 

8984 

854  • 

°27 1 

850  • 

°229 
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the  discharge  pulse  was  recorded  by  photographing  the  detector  waveform 
displayed  upon  an  oscilloscope  and  terminated  by  a 1 Kf2  resistor. 

D.  RESULTS 

The  small  signal  gain  go  is  defied  by  the  expccssion,  I/I  = e*p  (g  ( 
where  1 is  the  active  length  of  the  discharge  and  Io  and  I are  the  probe  las^r 
b.am  intensities  prior  to  and  after  propagation  through  the  discharge.  Gain 
profiles  as  a function  of  the  lower  rotational  level  P(J|  for  various  vibra- 
tional hands  are  shown  in  Figure  2.  Maximum  gain  of  approximately 
4 1/2  %/cm  was  measured  on  the  cp(10)  and  SPflO)  transitions.  Confirma- 
tion of  this  value  was  obtained  by  noticing  that  threshold  lasing  occurred 
by  reflection  off  the  cavity  windows  when  they  were  normal  to  the  optical 
axis.  Calculations  showed  that  the  cavity  loss  of  this  configuration  was 
approximately  4 1/2  %/cm.  At  threshold,  this  value  is  by  definition  equal 

to  the  gain.  During  the  gain  measurements,  the  windows  were  tilted  to 
prevent  oscillation. 

Most  of  the  data  points  within  a particular  vibrational  band  showed 
a smoothly  varying  peaked  distribution  as  shown  in  Figure  2 . The  5P(12), 

5 P(  1 3 ),  6P(10)  and  7P(10)  transitions  did  not  conform  to  this  pattern  and 
exhibited  anomalously  low  gain.  The  displacements  of  the  6P(10)  and 

7P(10)  data  points  from  their  corresponding  smooth  curves  are  indicated 
on  Figure  2 by  the  vertical  arrows. 

Lacina  and  McAllister (Z  1 have  attributed  the  absence  of  certain 
lasing  transitions  to  resonance  absorption  of  these  lines  by  overlapping 
transitions.  The  extent  of  the  overlap  anci  the  population  distribution 
determines  the  degree  of  absorption.  The  6P(10)  transition  is  overlapped 
by  the  8R(3)  line.  The  separation  of  the  unshifted  lines  is  0.  1 1 cm'1, 
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Figure  2 Measurements  of  Maximum  Small  Signal  Gain  as  Fui  cHons 

of  Terminal  Rotational  Level  for  Various  P Branch  Vibrational 
T ransitions 


and  for  the  present  experimental  conditions  the  lives  are  dominated  by 
pressure  broadening  with  a half-width  0.  03  cm'1(FWHM).  Thus,  absorp- 
tion occurs  due  to  overlap  between  the  wing  of  the  8R(3)  ransition  and  line 
center  of  the  f P(10)  line.  Additional  overlaps  exist  with  other  levels^2' 
which  are  not  expected  to  be  appreciably  populated.  The  peak  of  the  rota- 
tional distribution  is  expected  to  occur  close  to  J = 7;  therefore,  absorption 
will  mainly  occur  into  lower  rotational  levels.  The  coincidence  for  the 
5 P(  1 0 ) transition  which  occurs  with  the  SR(7)  line  is  even  closer,  to  within 
approximately  0.  03  cm’  \ This  very  close  coincidence  is  reflected  in  the 
absence  of  guin.  The  line  exhibits  rapid  total  absorption  of  the  probe  signal. 

Additional  support  for  the  hypothesis  of  resonance  absorption  by 
overlapping  transitions  was  derived  from  small  signal  gain  measurements 
on  the  5 P ( 1 3 ) transition  as  a function  of  E/N.  Figure  3 shows  four  oscillo- 
grams of  the  temporal  gain  behaviour  for  the  5P(13)  line  for  different 
values  of  E/N.  Notice  that  for  the  highest  E/N  the  line  exhibits  rapid 
total  absorption  and  that  as  the  E/N  is  reduced  the  absorption  becomes  less 
rapid  and  is  preceded  by  gain.  Eventually  the  probe  signal  is  incompletely 
absorbed.  The  electron  excitation  cross  sections  shift  towards  higher 
energies  for  higher  lying  vibrational  levels.  Consequently  reducing  E/N 
decreases  the  population  of  the  higher  levels  proportionality  more  than 
the  lower  levels.  Calculations  show  that  energy  transfer  by  V-V  collisions 
as  E/N  is  reduced  produces  similar  consequences.  Therefore,  as  E/N  is 
reduced  these  mechanisms  deplete  the  population  of  the  8R(7)  level 
correspondingly  more  than  the  5P(  1 3)  level  which  reduces  the  extent  of 
the  resonance  absorption  as  shown  in  Figure  3. 
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Detai's  of  the  modelling  calculations  have  been  given  elsewhere.'1’61 
Figure  4 shows  the  predicted  dependence  of  small  signal  gain  on  vibrational 
band  as  a function  of  P(J)  for  our  estimated  experimental  conditions.  In 
contrast  to  the  measurements  maximum  gain  is  predicted  to  occur  on  the 
1 2P(  1 0)  transition  and  the  magnitudes  of  the  peak  gain  on  all  transitions 
are  calculated  to  be  somewhat  higher  than  the  measured  values. 

Gain  saturation  by  the  probe  laser  beam  was  considered  but  re- 
jected as  a mechanism  offering  a possible  explanation  for  the  disagreement 
between  the  measured  and  calculated  values  of  maximum  gain.  Power 
output  measurements  on  the  20  liter  laser  set  a lower  limit  for  the  satu- 
rated flux  intensity  of  approximately  1 04  watcs/cm2.  After  amplification 
by  the  medium,  the  probe  beam  intensity  falls  three  orders  of  magnitude 
below  this  value.  The  lack  of  agreement  probably  reflects  uncertainties 
in  estimating  the  discharge  and  medium  properties  which  are  known  from 
previous  current  density  and  interferometric  measurements  to  show  con- 
siderable lack  of  uniformity.  ^1 

Resonance  absorption  of  the  6P(10)  and  7P(10)  transitions  by  the 
8R(3)  and  qR(3)  levels  respectively  was  included  in  the  gain  calculations. 

The  effect  is  indicated  in  Figure  4 by  the  vertical  arrows  which  indicate 
the  extent  of  the  gain  suppression.  The  predictions  are  in  good  agreement 
with  the  measurements  shown  in  Figure  2. 

Figure  5 shows  predictions  of  the  temporal  variation  of  small 
signal  gain  on  the  PM2)  transition  for  various  vibrational  bands.  The  dis- 
charge parameters  are  chosen  appropriately  to  describe  the  estimated 
experimental  conditions  and  are  indicated  on  the  figure.  Each  transition 
shows  absorption  initially  followed  by  gain.  This  is  consistent  with  the 
physical  model  envisaged  of  energy  being  sequentially  transferred  from 


IP*,  o/ 


level-to-level  up  the  vibrational  ladder  by  V-V  relaxation  collisions (b “ J 0) 
which  at  early  times  establish  absorption  and  at  later  times  produce  partial 
inversions  and  hence  gain.  The  small  gain  peak  predicted  initially  for  the 
6P(12)  transition  is  associated  with  direct  electron  pumping  of  the  v=6  level. 
Direct  excitation  is  far  less  significant  for  the  higher  lying  levels. 

For  comparison  with  the  calculations,  two  measurements  of  the 
temporal  gain  behavior  together  with  the  discharge  current  waveform  are 
shown  in  Figure  6.  The  9P(9)  oscillogram  is  representative  of  all  except 
the  resonantly  absorbed  levels  which  are  characterized  by  the  behavior 
shown  in  the  6 P(1  0)  oscillogram.  The  majority  of  transitions  do  not 
exhibit  the  initial  absorption  predicted  and  shown  in  Figure  5.  (The 
anomalous  transitions  exhibit  initial  absorption  followed  by  either  anoma- 
lously low  gain  or  further  absorption  at  later  times  in  contrast  to  the 
majority  of.  transitions  which  exhibit  gain  initially  which  decays  at  later 
times  to  zero  or  occasionally  to  absorption.  ) The  characteristic  trends 
of  the  temporal  gain  behavior  exhibited  in  Figure  5 such  as  the  times  to 
reach  threshold  and  peak  gain  as  a function  of  vibrational  level  were 
qualitatively  confirmed.  Figure  7 shows  the  experimental  and  predicted 
times  to  reach  threshold  gain  as  a function  of  upper  vibrational  level  for 
the  P(12)  transition.  Higher  lying  vibrational  levels  require  longer  times 
to  reach  threshold  gain  lending  additional  support  to  the  energy  transfer 
mechanism  proposed  for  establishing  partial  inversions. 

In  view  of  the  discrepancy  between  theory  and  experiment  regarding 
initial  absorption,  measurements  were  performed  which  confirmed  that  the 
detector  response  was  sufficiently  fast  to  respond  to  the  predicted  initial 
abso’-otion.  I urthermor e,  the  signal-to-noise  sensitivity  was  sufficient 

to  reveal  absorption  levels  as  low  as  0.  2%/cm,  well  below  the  predicted 
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Figure  6 


Temporal  Variation  of  the  Probe  Laser  Intensity  while 
Oscillating  on  the  6P  (10)  and  9P  (9)  Transitions.  The 
discharge  current  pulse  is  shown  at  the  center  of  the  figure. 
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values.  Using  the  predicted  population  distributions,  calculations  also 
confirmed  that  saturated  absorption  by  the  medium  was  not  responsible  for 
the  discrepancy.  Additionally,  further  kinetics  calculations  were  performed 
substituting  variously  increased  relaxation  rates  for  both  the  near  and  non- 
resonant collisions.  This  procedure  reduced  the  amplitude  of  the  predicted 
absorption  but  never  suppressed  it  entirely.  However,  there  is  no  justifi- 
cation for  this  procedure  since  be  rates  originally  employed^10*  agree 
well  with  the  limited  available  experimental  data.  (13)  Thus;  this  disagree- 
ment between  measurement  and  calculation  with  regard  to  initial  absorption 
is  not  understood. 

The  previous  measurements  were  performed  along  the  center  line 
of  the  discharge  volume.  In  order  to  investigate  the  effect  of  inhomogene- 
ities in  the  electrical  discharge  and  the  gas  density  the  probe  laser  beam 
was  translated  10  cm  to  the  edge  of  the  anode-cathode  region  while  remain- 
ing midway  between  them.  Similar  gain  measurements  showed  that  the 
character  of  the  temporal  variation  remained  generally  unchanged  but  did 
reveal  reduced  gain  of  approximately  25%  on  those  lines  investigated. 

In  conclusion,  the  measurements  show  good  qualitative  agreement 
with  the  kinetics  calculations  except  for  the  unexplained  absence  of  absorp- 
tion preceding  gain.  Differences  between  the  measured  and  calculated 
magnitudes  and  location  of  peak  gain  may  be  accounted  for  by  uncertainties 
associated  with  inhomogeneities  in  the  electrical  discharge  and  gas  density. 
The  measurements  support  the  suggestion  that  resonance  self- absorption 
exists  on  a number  of  transitions  and  calculations  which  include  this  effect 
compare  well  with  the  experimental  data.  It  is  recognized  that  this  mech- 
anism will  play  an  important  role  in  determining  the  atmospheric  trans- 
mission properties  of  the  high  pressure  carbon  monoxide  electric  discharge 


laser. 
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III.  UNSTABLE  RESONATOR 


An  ustable  resonator  configuration  has  been  designed  and  constructed 
for  use  in  conjunction  with  itie  water  vapor  absorption  cell  for  selected  line 
operation  of  the  20  liter  laser.  The  large  uniformly  filled  mode  volume 
offered  by  the  unstable  resonator  makes  it  a suitable  choice  for  such  an 
application. 

Power  measurements  and  burn  patterns  have  been  made.  The 
latter  indicates  a nonaniform  intensity  distribution  which  is  ascribed  to 
medium  inhomogeneities.  The  design  consisted  of  a single-ended  confocal 
configuration  employing  a square  convex  output  coupling  mirror  and  a 
larger  diameter  concave  mirror.  The  square  output  coupling  mirror 
offers  a more  uniform  intensity  distribution  in  the  lowest-loss  mode  and 
better  mode  discrimination  than  a circular  mirror. 

The  approximate  si zr>  of  the  square  output  mirror  was  determined 
by  the  area  of  the  laser  medium  which  was  to  be  utilized  (approximately 
5x5  cm)  divided  by  the  magnification  of  the  system.  The  exact  size  of  the 
output  coupling  mirror  was  arrived  at  by  requiring  that  the  equivalent 
Fresnel  number  of  the  system  be  half  integral.  The  equivalent  Fresnel 
number  F^  is  defined  by 

F = — — / 1 - — \ 

E 4 A L l 1 R J 

where  a is  the  radius  and  R the  radius  of  curvature  of  the  mirror,  L the 
mirror  separation  and  A the  wavelength.  The  purpose  of  this  constraint 
is  to  optimize  transverse  mode  discrimination. 
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It  is  generally  desirable  to  make  the  magnification  of  such  a system 
as  large  as  possible  in  order  to  minimize  the  obscuration  due  to  the  output 
coupling  mirror.  However,  if  it  is  made  too  large,  there  will  be  insuffi- 
cient feedback  to  sustain  lasing  action.  For  this  application,  ihe  coupling 
was  chosen  to  be  50^  which  corresponds  to  a magnification  of  \Z.  These 
considerations  determined  the  size  of  the  output  coupling  mirror  to  be 
3,  3 cm  square.  The  separation  of  the  mirrors  determined  by  the  size  of 
the  laser  and  the  water  vapor  cell  was  M-l/2  feet.  The  radius  of  the 
output  coupling  mirror  was  75  feet  and  the  larger  concave  mirror  100  feet. 
The  tolerances  imposed  on  the  figures  of  the  mirrors  were  chosen  on  the 
basis  of  the  "small  reduction  of  center  field  intensity  criterion.  " This 
required  that  the  figures  be  accurate  to  yg-.  The  two  elements  were 
arranged  on  the  center  line  of  the  discharge  and  carefully  aligned  with  an 
autocollimato  r. 

A photograph  of  the  near  field  burn  pattern  obtained  is  shown  in 
Figure  8,  the  location  of  the  anode  and  cathode  are  indicated.  Also  visible 
in  the  photograph  is  the  outline  of  the  square  output  coupling  mirror  and  in 
the  lower  portion  one  of  the  arms  of  the  spider  mount  support.  The  highest 
intensity  occurs  near  the  lower  portion  of  the  pattern  which  is  closest  to 
the  cathode.  Similar  measurements  by  Northrop  yielded  quite  different 
results  with  areas  of  highest  intensity  being  nearer  to  the  anode.  The 
lack  of  uniformity  in  their  case  was  interpreted  in  terms  of  electron 
scattering  which  reduced  the  relative  electron  density  near  the  anode 
thereby  increasing  the  medium  resistivity  and  hence  the  energy  deposition. 
In  the  present  system  this  effect  is  also  expected  to  occur  but  the  intensity 
distribution  is  believed  to  be  dominated  by  the  density  variation  within  the 
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Figure  8 The  Near  Field  Burn  Pattern  Obtained  Using  the  Unstable 
Resonator  Configuration 
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medium.  Previous  interferometric  measurements  of  medium  homogeneity 
revealed  considerable  density  gradients.  The  cryogenic  gas  entering 
through  the  porous  cathode  is  heated  through  conduction  and  radiation  as 
it  traverses  the  cathode  - anode  gap.  However,  since  ts3  anode  is  cooled 
by  circulating  LN^  the  temperature  of  the  gas  in  the  vicinity  of  the  anode 
is  reduced  which  imposes  a reverse  temperature  gradient.  The  net 
effect  is  to  relatively  reduce  the  E/N  in  the  upper  discharge  region  towards 
the  anode  thus  reducing  the  energy  deposition  and  producing  the  nonuniform 
pattern  shown  in  the  figure.  The  lateral  displacement  of  the  pattern  is 
believed  to  be  associated  with  misalignment  of  the  unstable  resonator 
optics. 


IV.  WATER  VAPOR  ABSORPTION  CELL 


The  water  vapor  absorption  cell  is  to  be  used  in  conjunction  with 
the  unstable  resonator  in  an  intracavity  configuration.  The  purpose  is  to 
pursue  studies  of  selected-line  operation.  A schematic  of  the  absorption 
cell  and  the  associated  equipment  is  shown  in  Figure  9.  The  construction 
is  of  stainless  steel  except  for  the  oven  enclosure  which  is  thin-walled 
aluminum. 

The  oven  enclosure,  transfer  line  and  water  reservoir  are  wrapped 
with  heater  tapes,  with  a number  of  independent  temperature  controllers. 
Temperature  monitoring  is  provided  by  ten  thermocouples  placed  at  vari- 
ous points  over  the  system.  The  pressure  within  the  cell  is  measured 
with  a high  temperature  differential  pressure  transducer.  Additive  gases 
may  be  introduced  to  the  cell  via  the  metering  valve.  The  assembly  is 
pumped  by  a 3 cfm  mechanical  pump  which  permits  a base  pressure  of 
approximately  10yi  to  be  achieved.  A liquid  nitrogen  cold  trap  permits 
the  water  vapor  to  be  trapped  whenever  it  is  necessary  to  evacuate  the 
system.  The  clear  aperture  cf  the  cell  is  5 l/2"  in  diameter,  it  is  limited 
to  3 l/2"  by  the  two  cylinders  at  either  end  of  the  oven  enclosure.  The 
purpose  of  the  cylinders  is  to  reduce  convective  air  circulation  around  the 
windows.  The  optical  length  of  the  cell  is  20".  The  cell  windows  are  of 
calcium  fluoride  and  6"  in  diameter  they  are  placed  at  normal  incidence  to 
the  optical  axis  and  have  been  anti- reflecting  coated  at  5 M- 
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Figure  9 Schematic  of  the  Water  Vapor  Absorption  Cell 
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Prior  to  installing  the  cell  inside  the  unstable  resonator  cavity 
interferometric  measurements  will  be  performed  to  check  the  medium 
quality  inside  the  absorption  cell. 


V.  NEW  DISCHARGE  CAVITY/FLOW  SYSTEM 


A new  transverse  gas  uow  system  has  been  designed  and  constructed 
in  this  continuing  program  to  overcome  the  thermal  nonuniformities  ob- 
served in  the  original  discharge  cavity.  A schematic  on  the  flow  system 
is  shown  in  Figure  10.  The  basic  concept  of  the  new  system  is  to  intro- 
duce the  gas  through  a heat  exchanger  and  allow  free  convection  boundary 
layers  to  be  established  on  the  room  temperature  side  and  end  walls.  The 
flow  direction  is  vertically  upwards  to  ensure  that  buoyancy  effects  are 
favorable  and  to  avoid  recirculation  zones.  The  transverse  flow  configura- 
tion in  this  new  cavity  design  has  the  advantage  of  utilizing  the  smallest 
cross  sectional  area  in  the  system  resulting  in  maximum  flow  velocity  for 
a given  pumping  speed. 

The  new  cavity  flow  system  is  designed  for  operation  over  the 
temperature  range  from  80°K  to  300°K  and  the  pressure  up  to  1 atm.  The 
overall  design  goal  is  to  obtain  the  beam  quality  that  is  1.  5 times 
diffraction-limited  or  better.  This  corresponds  to  a maximum  of  a 1 / 1 0 
fringe  variation  across  the  cavity.  A new  flow  system  is  designed  to  be 
compatible  with  the  existing  20  cm  x 100  cm  electron  gun  and  to  operate 
with  a maximum  discharge  width  of  the  order  of  10  cm.  The  system  is 
designed  to  operate  initially  with  a transverse  flow  velocity  of  100  cm  per 
second  and  a corresponding  mass  flow  cf  the  order  of  150  grams  per 
second.  Provision  has  been  made  for  increasing  the  transverse  flow 
velocity  up  to  1000  cm  per  sec  should  this  prove  necessary. 
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Figure  10  Schematic  of  Transverse  Flov'/Discharge  Cavity  System 
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The  design  concept  of  the  heat  exchanger  was  to  use  a pre-rooled 
porous  plug  with  sufficient  thermal  capacity  to  cool  the  gas  for  transient 
operation  of  approximately  10  sec.  This  is  more  than  sufficient  for  pulsed 
operation  and  avoids  the  necessity  for  construction  of  a massive  heat 
exchanger  with  all  its  associated  liquid  nitrogen  plumbing. 

A.  BOUNDARY  LAYER  CALCULATIONS 

Medium  homogeneity  requirements  for  a diffraction- limited  output 
beam  place  stringent  limits  on  the  allowable  variations  in  optical  path 
length  in  the  cavity.  Since  the  cavity  walls  are  warmer  than  the  gas, 
buoyancy  effects  appear  along  with  forced  convection.  In  addition,  a uni- 
form flow  in  velocity  entering  the  cavity  is  required  to  allow  the  formation 
of  boundary  layers  on  the  windo.vs.  Boundary  layers  (as  opposed  to  re- 
circulation zones)  are  desirable  because  they  concentrate  the  density 
change  into  a narrow  distance,  leading  to  a low  path  length  variation  rela- 
tive to  a uniform  medium.  A configuration  such  as  shown  in  Figure  10  is 

clearly  suitable  for  this  purpose. 

The  maximum  variation  in  the  quantity  AjC  , where 

A P " Poo  , m 

AC  = / — dy  (1) 

o 

between  any  two  paths  will  be  set  here  at  l/lO  of  a fringe  corresponding  to 
a l.  5 times  diffraction  limited  central  spot.  In  the  above,  p is  local 
density,  is  the  bulk  flow  mean  density,  I is  the  cavity  length,  and  y is 
distance  perpendicular  to  the  windows  in  the  lasing  direction.  At  with 

l/2  amagat  density,  the  value  of  A£  for  l/lO  fringe  is  3 mm.  The  boundary 
layer  calculations  given  below,  indicate  that  this  constraint  can  be  readily 
satisfied. 
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There  are  three  dimensionless  parameters  of  interest.  The  first 

one  is  the  Reynolds  number,  Red  = ud/v,  where  u is  the  bulk  flow  velocity, 

d is  he  cavity  width,  and  v is  the  kinematic  viscosity  of  the  gas.  The 

Reynolds  number  is  the  parameter  associated  with  forced  convection.  For 

u=  1°Z  cm/sec  (a  lower  bound)  d=  10  cm,  v = 0.  2 cm2/sec.  Re  = 5 x 1 03 

d 

Free  (buoyancy-driven)  convection  is  controlled  by  the  Grashof  number 
Grd  = gP  0 wd  /v  , where  g is  the  acceleration  of  gravity,  (3  is  the  com- 
pressibility of  the  gas  and  is  the  temperature  difference  between  the 
window  and  the  bulk  flow.  With  Q ^ = 2 00°K  (bulk  flow  at  100°K),  and  an 
effective  P = ^ where  T*  = 2 00°K,  Grrf  = 2.  5 x 1 07. 

The  parameter  that  controls  the  relative  importance  of  forced 
convection  and  buoyancy  is  the  Froude  number  Frd  = u/(g|3  Q d)1/2.  The 
Froude  number  is  simply  the  ratio  Re^/Gr  ,l/2.  In  our  case,  Fr,  =*  1 
Buoyancy  is,  therefore,  not  a small  perturbation  on  forced  convention,  but 
rather  of  equal  importance.  Because  of  this,  it  is  essential  to  use  the 

optimum  flow  configuration  of  Figure  1 for  which  the  buoyancy  forces  are 
in  the  flow  direction. 

The  Reynolds  and  Grashof  number  can  be  written  in  terms  of  dis- 
tance along  the  wall  (x)  to  determine  whether  the  boundary  layer  will  be 
laminar  or  turbulent,  at  x = 40  cm,  Rex  = 2 x 1 04  and  Gr  =1.6x1  09. 

A pure  forced  laminar  boundary  layer  will  undergo  transition  to  turbulence 

5 

at  Rex  10  , whereas  a free  layer  will  do  so  at  Grrf  ~ 1 09.  There  is  a 
possibility,  then,  that  the  boundary  layer  on  the  window  may  be  turbulent. 

The  change  in  optical  path  length  A £ in  the  boundary  layer  is  ob- 
tained by  integration  of  Eq.  (1)  from  0 to  6,  where  6 is  the  boundary  layer 
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thickness.  The  change  in  path  length  which  is  caused  by  turbulent  fluctu- 
ations in  the  boundary  layer  is  much  smaller  than  A£  in  Eq.  (1),  and  will 
be  ignored.  If  A£  did  not  vary  from  point  to  point  (in  the  stream-wise 
direction)  in  the  boundary  layer,  each  ray  would  travel  the  same  path 
length.  It  is  the  change  in  Af  from  point  to  point  as  the  boundary  layer 
becomes  thicker  that  affects  the  medium  homogeneity. 

The  problem,  then  is  the  computation  A£  as  a function  of  stream- 
wise  location  for  a boundary  layer  on  a vertical  heated  wall  where  both 
forced  and  free  convection  are  present.  The  structure  of  the  laminar  layer 
has  been  obtained  analytically  and  numerically  by  several  authors.  ^ ^ 

In  this  report  an  approximate  solution  is  obtained  for  the  first  order  per- 
turbation on  the  asymptotic  solutions  for  pure  forced  convection  and  pure 
free  convection  using  an  integral  formulation  for  the  turbulent  boundary 
layer  equations.  Characteristic  length  scales  for  optical  path  variations 
will  be  presented. 

1.  Forced  Convection  Asymptotic  Case 

In  the  absence  of  pressure  gradients  and  body  forces  the  boundary 

. (19) 

layer  momentum  equation  for  the  pure  forced  convection  case  is 


J u K ' u 


) dy  = T 


w 


(2) 


where  6j  is  the  pure  forced  convection  boundary  layer  thickness,  u is 
velocity,  u^  is  velocity  far  away  from  the  wall,  and  7^  is  wall  shear 
stress. 
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As  suming 


• ■ •-  (4 


1/7 


Tw  = 0.  0228  p 


1/4 


(3) 


the  well  known  result  for  5j  is  obtained, 


= 0.  376 


(4 


U\^  4/5  _ 0.  376x 

’ (R'x>1/5 


(4) 


where  v is  the  kinematic  viscosity,  x i,s  distance  along  the  window,  and 
Ke  is  Reynolds  number,  Re  = u x/v . 

If  the  boundary  layer  is  heated  in  a gravity  field,  buoyancy  forces 
arise.  If  the  buoyancy  force  is  included  in  the  momentum  equation  using 
the  approximation  A p/p  = (T  Tj/T^,  Eq.  (2)  be< 


icomes 


5i  +<Y 


/ u (u_  - u)  dy  = Tw  - f 


51  * V 


9 dy 


(5) 


where  g is  the  acceleration  of  gravity,  is  the  gas  temperature  far  away 
from  the  wall,  0 = T - T^,  and  6/  is  the  perturbation  in  boundary  layer 
thickness  caused  by  the  addition  of  buoyancy. 


9 = 9»  [‘  ’ («)  /7]  <6) 
(where  ©w  is  - T^).  We  will  also  assume  that  the  buoyancy  force  is 

small,  and  that  6/  « 6j.  It  is  possible  to  derive  equations  for  parameters 

describing  the  perturbations  in  boundary  layer  velocity  and  temperature 

from  the  momentum  and  energy  equations  taking  constant  6^.  This 
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procedure  requires  that  one  choose  suitable  profiles  for  the  perturbations. 
Since  there  is  no  information  available  on  such  profiles,  any  profile 
choice  would  necessarily  be  arbitrary.  For  this  reason,  we  will  assume 
that  buoyancy  appears  as  a change  5j'  in  the  boundary  layer  thickness 
with  the  unperturbed  velocity  profile  (Eq.  (3)).  Then  we  can  obtain  an 
equation  for  the  perturbation  in  caused  by  buoyancy  by  use  of  Eqs. 

(2)  - (6)  and  this  leads  to  the  result 


5 


1 


2 

u 

oo 


(7) 


Buoyancy  thus  reduces  the  boundary  layer  thickness  at  a given 
value  of  x,  since  the  perturbed  thickness  is  ( 6 j + ^j'). 

2.  Free-Convection  Asymptotic  Case 

Eckert  and  Jackson^^  have  solved  the  problem  of  pure  free- 
convection  turbulent  boundary  layer.  Taking  integral  equations  for  mo- 
mentum and  energy,  assuming  Pr  = 1,  and  using 


one  can  obtain 

52  = 0.  5b5  (Grx)_1/10  x (10) 

where  is  the  free  convection  layer  thickness  and  Gr  is  the  Grashof 

3 2 

number,  Gr^  = g (G^/T^)  x /v  . For  comparison,  in  the  laminar  case 
£>2  ~ x1/4  instead  of  ~ x‘ ' 1 ^ as  in  Eq.  (10). 
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Introducing  a perturbation  consisting  of  an  external  flow  (velocity  u ) 

OO 

into  the  integral  energy  equation,  we  obtain 


6Z  + 62' 


pCPs  / 

0 


(u  + u ) 0 dy  = q 


w 


(11) 


Proceeding  as  before,  we  obtain  an  equation  for  the  perturbation, 
62\  in  62  caused  by  the  external  flow.  This  leads  to  the  solution 

u /g  (0  /T  )\'6/10 

«2'  * -0^'0-f  ( V2-)  *2/l°  (12) 


The  external  flow  reduce  6^  at  a given  x,  since  the  perturbed 
thickness  is  (§2  + ' ). 

3.  Results  for  Optical  Path 

The  results  obtained  above  for  the  perturbed  boundary  layer  thick- 
ness can  now  be  used  to  calculate  optical  path  lengths  through  the  boundary 
layer. 

With  the  approximation  (p  - p ^ )/p^  = a/T^  , Eq.  (1)  becomes 

6 

f 9 , 

= I — ay  (13) 

0 00 


Af 


Using  Eqs.  (4),  (6),  (7)  and  (13),  we  obtain  for  the  forced- 
convection  asymptotic  case 

/ 

Gr 


Af 


L w 0.  04  7 ~ x Re  ’l/5  [1-0. 

00  x V 


64 


Re 


(14a) 


Similarly,  using  Eqs.  (9),  (10)  and  (12) 
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if  i \J 


(15a) 


Af,  « 0.  074  ■==  x Gr 
2 T x 

CO 


{ TT) 


The  dominant  parameter  in  this  problem  is  a Froude  number  based 
on  x,  where 

u Re 


Fr 


.«  <ew/T»>  *]  1/2  ' lG’-)I>'/2 


(16) 


There  is  then  a scale  length  L#  for  distance  along  the  window, 

2 


L*  = 


g (9  /T  ) 
^ w 00 


(17) 


If  x « L*,  forced  convection  is  dominant.  As  x becomes  very  large 
(x  » L*),  free  convection  (i.  e.  , buoyancy)  dominates.  Neither  Eq.  (14) 
nor  Eq.  (15)  applies  when  x/L*  ~ 1. 

In  nondimensional  form,  the  above  results  are 


,4/5 


Af  * « 0.  047  x*  ' J (1  - 0.  64  x*) 


(Forced  Convection)  (14b) 
Af  * « 0.  074  x*'/,0G  !.53x*'1//Z)  (Free  Convection)  (15b) 


Here 


x*  = x/L* 


Af*  = AC/Y* 


TURB 

where  ^jjRB  *8  a sca^nS  length  for  the  optical  path  variation,  and  is 
defined  by 


(18) 

(19) 


(9  /T  )>/5  u 7/5  v’/5 

w 00  OO 
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(20) 


,1 

I 


Following  a procedure  similar  to  that  presented  above  for  a laminar 
layer  leads  to  a different  scaling  length 


YL\M 


(21) 


while  L*  is  as  in  Eq.  (17).  Values  for  L*,  and  YBAM  as  functions 

of  u and  0 /T  for  CO,  N_  mixtures  at  1 amagat  density  are  shown  in 

00  W OO  £ ‘ ' 

Figure  2.  Since  the  values  of  v for  CO  and  are  nearly  identical,  the 
mixture  ratio  is  unimportant.  In  Figure  11  (9w/Tm)  ?nd  v have  been  evalu- 
ated using  a reference  temperature  halfway  between  and  T to  replace 
T for  large  (T  - T ).  Regions  in  x,  u , and  T where  the  flow  is 

OC  W OO  00  OO 

5 q 

turbulent  have  been  delineated  using  the  criteria  Re  =2x10  or  Gr  = 10  . 

x x 

Using  the  results  of  Ref.  15,  A£  for  a laminar  layer  can  be  obtained. 
Figure  12  shows  a plot  of  Af  vs  downstream  distance  x from  the 
leading  edge  (the  heat  exchanger)  using  a reference  temperature  T*  of 
200°K  to  replace  T , for  conditions  of  interest  for  the  CO  laser, 
u^  = 100  cm/sec,  d = 10  cm,  and  0w  = 2 00°K,  for  a CO-N2  mixture. 

Figure  13  shows  A£  vs  x for  a turbulent  boundary  layer  for  the 
same  operating  conditions  as  above.  Since  the  asymptotic  expressions 
Eqs.  (14)  and  (15)  do  not  join  up  smoothly,  and  are  not  valid  beyond  the 
point  where  the  perturbations  are  — 25%  of  the  unperturbed  value,  lines 
were  drawn  tangent  to  the  perturbed  free  and  forced  Af  at  the  point  where 
the  perturbation  is  25%  of  the  unperturbed  valve.  These  lines  form  the 
combined  flow  line  shown  in  Figure  13.  Lines  showing  the  perturbed 
valves  (Eqs.  (14)  and  (15))  are  also  included  in  Figure  13. 
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A I (cm) 


From  Figure  12  it  can  be  seen  that  if  the  window  extends  from 

x = 20  cm  to  x = 40  cm;  and  the  flow  remains  laminar,  the  variation  in  Af 

from  top  to  bottom  will  be  0.  07  cm,  or  1.  4 mm  for  two  windows.  On  the 

other  hand,  assuming  that  the  boundary  layer  is  laminar  at  x = 20  cm 

(Rex  and  Gvx  are  well  below  transition  levels  at  x = 20  cm)  and  turbulent 

at  x = 40  cm,  from  Figures  12  and  13  the  variation  in  Af  would  not  be  any 

greater.  Since  the  boundary  layer  is  transitional,  there  is  little  difference 

in  the  thicknesses  which  is  well  w-ithin  the  required  limit.  The  turbulent 

layer  grows  at  a very  fast  pace,  however,  so  that  for  larger  systems,  say 
2 

x = 10  cm,  the  difference  between  a turbulent  layer  and  a laminar  layer 
would  be  large. 

B.  DESIGN  OF  HEAT  EXCHANGER 

Following  some  preliminary  investigations  it  was  decided  to  use 
porous  sintered  bronze  plates  for  the  heat  exchanger.  This  choice  was 
based  upon  the  requirement  of  a large  thermal  capacity  and  efficient  heat 
transfer  to  the  gas.  It  should  be  noted  that  free  stream  temperature  vari- 
ations due  to  dynamic  effects  are  negligible  in  the  present  case  with  flow 
velocities  less  than  1000  cm  per  sec. 

The  overall  design  of  the  heat  exchanger  is  shown  in  the  schematic 
of  Figure  14.  It  consists  basically  of  three  uniform  sintered  bronze  plates 
mounted  in  an  isothermal  aluminum  box.  A series  of  1 /4  inch  copper 
tubes,  spaced  2 inches  apart,  are  embedded  in  the  sintered  bronze  plates 
and  connected  via  manifolds  to  the  coolant  supply,  usually  liquid  nitrogen. 
The  manifolds  are  split  into  two  counterflow  sections  to  reduce  any  trans- 
verse temperature  gradients.  The  porosity  of  the  plates  was  chosen  to 
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provide  a pressure  drop  of  approximately  1 psi  across  each  plate  based 
on  laboratory  measurements. 

It  can  be  shown  analytically  that  large  scale  inhomogeneities  pres- 
ent in  the  laminar  flow  issuing  from  the  porous  plates  decay  very  slowly 
for  the  Reynolds  numbers  of  interest.  Thus,  it  is  necessary  to  insure 
that  the  flow  leaving  the  porous  plates  has  good  initial  temperature  uni- 
formity, a fraction  of  a degree,  to  prevent  significant  contribution  to  the 
quantity  A*  defined  in  Eq.  (1).  Calculations  were  made  to  estimate  the 
thermal  accommodation  of  the  gas  to  the  porous  plates  to  ensure  that  the 

gas  leaves  the  heat  exchanger  at  a uniform  temperature.  These  calcula- 
tions are  outlined  below. 

The  porous  plates  can  be  viewed  as  a collection  of  long,  narrow 
tubes  of  diameter  D.  The  heat  transfer  in  such  a tube  is 

n C u AT  - 4 a — ' x - 4 /v  ^ T \ Ax 

P p - 4Vd"  - 4^j  r (22) 

where  um  is  the  mean  flow  velocity  AT  is  the  temperature  chance,  Ac  is 

distance  along  the  tube,  qw  is  the  heat  transfer  rate  per  unit  area,  and  r is 
tube  radial  distance.  The  pressure  drop  has  the  form 

Ap  = 4 r - 4 9u  \ Ax 

w D y1  Qr  J D (2  3) 

where  7 w is  the  wall  shear  stress.  If  Pr  - 1,  then  the  T(r)/9  and 

w 

u(r )/um  should  be  roughly  the  same,  so  that  Eqs.  (22)  and  (23)  can  be 
combined  to  give 

^ 1 dp 

6 dx  ~ 2 dx  (24) 

w Pu™ 
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Since  the  pressure  drop  across  the  plug  can  be  ~ 1 psi  (~  1 0^  dynes/ 
cm  ) and  pu2  ~ 1 0 J x (100)2  ~ 10  dynes/cm2,  the  distance  for  the  gas  to 
accommodate  to  the  wall  temperature  would  be  ~(10/lo'’)  x (thickness  of 
the  plug). 

1.  Required  Plate  Thickness 

It  will  be  convenient  to  assume  that  the  plug  is  made  up  of  a number 
of  thin  plates  each  of  thickness  L.  It  will  also  be  assumed  that,  as  shown 
above,  the  gas  and  metal  in  each  plate  accommodate  to  each  other.  All 
the  plates  are  pre-cooled  to  the  desired  final  gas  temperature  before 
flowing  the  lasing  mixture.  The  flow  is  then  turned  on  for  a period  ~ 1 - ] 0 
sec.  To  ensure  that  the  gas  leaves  the  stack  with  the  required  uniformity, 
the  last  plate  must  not  rise  in  temperature  more  than  ihe  required  uni- 
formity level  during  the  run. 

Flowing  a gas  of  initial  temperature  T^  through  a thin  plate  (i.  e.  , of 
uniform  temperature)  initially  at  Tw  for  a time  At  = L/u,  with  a final 
temperature  of  T^.,  a heat  balance  can  be  written  lc  the  gas  and  the  plate, 

Pr  r <T  - Tf>  = pw  C (Tf  - Tw)  (25) 

w 

In  the  above,  p is  the  gas  density,  A is  the  plate  area,  C is  the  gas  heat 

capacity,  p^  is  the  plate  density,  and  C is  the  plate  heat  rapacity. 

pw 

Taking  rdT/dt)plate  = (Tf  - Tw)/(L/u),  and  pw  C » p C a differential 
equation  for  the  temperature  of  the  first  plate  T}  can  be  obtained  from 
Eq.  (2  5) 
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(27) 


where  Tg^,  is  the  temperature  of  the  gas  entering  the  stack  of  plates,  and 

Tw  is  the  initial  temperature  of  the  stack.  For  the  second  slab,  the  incom 
o 

ing  gas  temperature  is  no  longer  T , but  the  gas  temperature  leaving  the 
first  plate,  or  thus, 

dT2  (Tx  - T ) 

~dT  * r (28) 

and  so  on  for  T3 T^.  The  mass  flow  through  each  plate,  p u A 

(where  A is  the  plate  area)  is  the  same,  and  the  area  of  each  plate  is  the 
same,  so  that  T is  the  same  for  each  plate. 

It  can  be  shown  by  induction  that  the  temperature  of  the  n-th  plate 
T in  the  stack  with  time  is 


T - T 
n 8, 


n-1 


T - T . . 

w g ^ 

o 6o  m=0 


[t  ' (n'1}  jf]  - [t  - (n-1)  £]/t 

. _m  e 


m ! r 


(29) 


With  typical  values  (L  ~ 1 inch  = 2.  5 cm,  u = 100  cm/sec), 

T ~ (5  x 0.  4)  x 2.  5 cm/ (1  0 3 x 1)  x 1 02  cm/sec  ~ 50  sec.  Then  for  a 

10  - second  run,  with  5 plates  and  T = 3 00°K  T = 1 00°K 

’ w * 

o o 

= (100  + 0.  06  )°K 


(30) 


t = 1 0 sec 

Because  of  frictional  heating  in  the  plates,  the  absolute  temperature 
of  the  whole  stack  will  rise  during  the  run.  The  energy  deposited  in  the 
plates  is  A puAt.  This  corresponds  to  a temperature  rise  AT, 
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AT  ^ 


(31) 


For  Ap/L  ~ 1 psi/inch  ~2.  5 x 10^  dyne/cm^/cm,  AT  < 1.  2°K  over 
a 10  sec  run,  or  0.  12  I</sec.  If  the  flow  uniformity  is  5%,  then  the 
temperature  variation  will  he  0.  1 x 1.2  ^ l/8°K.  There  is  also  a differ- 
ence in  temperature  between  the  center  of  each  pore  and  the  outer 
part  ~ (y  - 1 ) PrM  /8  — 10  where  M is  the  flow  Mach  number  or 

AT  ~ 10  x 100°K  ~ 10  ^ K,  which  is  insignificant. 

In  summary,  if  the  plates  are  properly  precooled,  the  gas  will 
leave  the  stack  with  variations  in  temperature  less  than  the  level  of  l/6°K 
required  for  diffraction-limited  operation. 

C.  DISCHARGE  ELECTRODE  CONFIGURATION 

Some  consideration  has  been  given  to  th  *sign  to  the  discharge 
electrodes  in  order  to  minimize  the  nonuniformities  in  the  deposition  of 
electrical  energy  in  the  gas.  Th^se  nonuniformities  arise  from  scattering 
of  the  primary  electrons  and  nonuniform  transverse  electric  field  distribu- 
tion due  to  the  finite  extent  of  the  electrodes.  These  in  turn  give  rise  to 
varying  conductivity  in  the  discharge  region. 

The  discharge  can  be  thought  of  as  nonuniform,  nonlinear  two- 
dimensional  resistor  between  the  anode  and  cathode.  Mathematically  the 
resistor  is  described  by  Poisson's  equation 

V.  = 0 (32) 

E = -V<£  (33) 

where  <f>  is  the  potential  and  E is  the  electric  field.  Eq.  (32)  is  nonlinear 
because  the  conductivity  a is  a function  of  E. 
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(2  1 ) 

Equation  (32)  has  been  studied  analytically  for  special  cases. 

This  was  done  by  a perturbation  expansion.  The  small  parameter  in  the 
expansion  was  the  ratio  E /E^  where  E^  is  the  electric  field  parallel 
to  the  anode/cathode  surfaces  and  E^  was  the  electric  field  perpendicular 
to  the  anode/cathode  surfaces.  It  was  further  assumed  that  the  anode  and 
cathode  were  planar  and  infinite.  However,  the  E-beam  was  finite  in 
extent.  For  a practical  laser  it  was  necessary  to  solve  Eq.  (32)  in 
arbitrary  geometry  to  allow  for  changes  in  design.  Because  of  the  com- 
plexity of  the  problem  it  is  necessary  to  seek  a numerical  solution  to 
Eq.  (32).  This  was  done  by  modifying  JASON,  a nonlinear  electro- 
static code  developed  by  Dr.  Sackett  of  Lawrence  Berkeley  Laboratory. 
This  code  has  been  tested  to  insure  that  it  does  in  fact  solve  nonlinear 
problems  accurately  with  errors  in  the  potential  and  electric  field  of  less 
than  0.  1%  and  1%  respectively. 

The  conductivity  a is  a function  of  the  local  electric  field  and 
ionization  rate  S by  the  high-energy  electrons.  This  dependence  will  vary 
from  one  gas  mixture  to  the  next.  No  measurements  of  the  recombination 


rate  have  been  made  in  an  80-20  N7-CO  mix.  However,  measurements 


(23) 


in  pure  N^  indicate  that  the  recombination  rate  varies  as  (Tg) 


-1.  8 


A similar  dependence  on  the  electron  temperature  Tg  was  obtained  for 


CO. 


(24) 


In  the  subsequent  analysis  it  is  assumed  that  the  recombination 

- 1 8 

rate  in  the  mixture  of  interest  is  proportional  to  (Tg)  ’ . From  the 

AERL  Boltzman  code  it  is  known  that  T in  an  80-20,  N--CO  mix  is 

e 2 

proportional  to  (E)^’  ^ while  the  drift  velocity  v^  oc  (E)^’ Combining 

(8 ) 

this  information  one  can  write 
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n - C(E)‘0,2S0,5 


where  C is  a constant. 


The  ionization  rate  S is  proportional  to  the  energy  deposition  by 
high-energy  electrons.  For  the  conditions  in  the  CO  laser  this  can  be 


written  approximately  as 


S ~ c x n 


1 fan-1  (lit)  - .an-1  (*^±) 


where  y = cE/2w  and  w is  the  electron  energy.  The  electron  beam  is 
assumed  to  have  a width  2a  located  at  x = 0 and  between  -a  < v :<  a 

Figure  15  shows  the  computer  plot  for  lines  of  constant  potential. 
Note  that  the  lines  are  more  closely  spaced  at  the  edge  of  both  electrodes. 
In  these  regions  the  electric  field  is  about  25%  larger  than  the  field  in  the 
center  of  the  laser  cavity.  Figure  16  illustrates  the  result  of  assuming 
that  the  ionization  rate  is  given  by  Eq.  (35).  Notice  that  the  lines  of 
constant  potential  are  not  evenly  spaced  in  Figure  16.  Con. paring 
Figures  15  and  16  it  is  evident  that  while  near  the  anode  the  lines  in 
Figure  16  are  more  closely  spaced,  the  reverse  is  true  near  the  cathode. 
The  electric  field  at  the  anode  Is  30-40%  greater  than  at  the  cathode.  This 
results  in  a 30-35%  increase  in  the  electrical  energy  deposited  near  the 
anode  as  shown  in  Figure  16.  The  highest  rate  of  energy  deposited  in 
Figure  17  is  in  the  region  where  the  anode  surface  begins  to  curve.  This 
is  a result  of  a slight  bulge  in  the  anode  boundary.  Such  a bulge  is  due  to 
the  method  the  computer  maps  the  boundary  and  is  not  present  ’ ^ the  actual 
cavity.  Finally,  the  optical  gain  is  expected  to  have  the  same  topology  as 
the  electrical  energy  deposited. 
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UNNt 


D.  FABRICATION  AND  TESTING 


r ' 


The  entire  flow  system  and  heat  exchanger  unit  have  been  con- 
structed, assembled  and  integrated.  Integration  with  the  electron  gun  will 
not  take  place  until  after  the  flow  and  homogeneity  tests.  The  electrode 
system  will  be  fabricated  in  the  coming  month  and  will  be  installed  after 
the  initial  flow  tests. 

Before  the  final  assembly  of  the  heat  exchanger  a series  of  thermo- 
couples was  embedded  in  the  central  porous  plate  in  order  to  measure  the 
temperature  distribution  during  the  cooling  process.  Heating  coils  were 
mounted  underneath  the  heat  exchanger  to  allow  for  rapid  (several  hours) 
warmup  of  the  system.  The  net  mass  of  the  heat  exchanger  is  the  order 
of  105  grams  and  a heat  input  of  the  order  of  1 0 joules  is  required  to 
regain  room  temperature. 

Preliminary  cooling  tests  with  liquid  nitrogen  resulted  in  the 
development  of  several  leaks  in  the  liquid  nitrogen  manifold  system.  This 
was  presumably  due  to  the  very  rapid  contraction  and  resulting  mechanical 
strains  as  liquid  nitrogen  entered  the  system.  To  avoid  this,  a cooling 
cycle  has  been  developed  to  allow  the  fairly  uniform  and  initially  much 
slower  cooling  of  the  system.  Tests  were  then  conducted  with  the  entire 
system  being  cooled  down  to  liquid  nitrogen  temperature.  Thermocouple 
measurements  indicated  the  maximum  variation  of  2°C  at  liquid  nitrogen 
temperature  across  the  entire  central  plate  of  the  heat  exchanger. 

There  is  still  some  evidence  of  leaks  developing  at  or  near  liquid 
nitrogen  temperature  in  the  manifold.  These  leaks  only  develop  at  very 
low  temperature  and  are  not  apparent  at  room  temperature.  They 
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presumably  occur  due  to  mechanical  stress  in  the  Swagelock  fittings  used 
to  construct  the  liquid  nitrogen  manifold.  Although  these  leaks  will  not  hold 
up  flow  testing  in  the  immediate  future,  provision  is  being  made  for 

replacement  of  the  Swagelock  fittings  with  an  all  welded  manifold  to  try 
and  avoid  leaks  in  the  future. 
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